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Abstract 



Oh | It is shown that the high-energy expansion of the scattering ampli- 

jy^ • tude calculated from Feynman diagrams factorizes in such a way that 

. it can be reduced to the eikonalized form up to the terms of inverse 

\ power in energy in accordance with results obtained by solving the 

Klein-Gordon equation. Therefore the two approaches when applied 
to the suppression of the emission of soft photons by fast charged par- 
00 | tides in dense matter should give rise to the same results. A particular 

limit of thin targets is briefly discussed. 

o ' 

1 Introduction 

5 : 

i-O ' In 1953, Ter-Mikaelian [p]] noticed that at high energies the longitudinal mo- 

mentum transferred at each scattering to an electron traversing a medium 
becomes very small. It enlarges the effective formation length for emitted 
^ ■ photons and leads to specific effects in crystals. Soon, Landau and Pomer- 

>• . anchuk remarked that it would result in the suppression of radiation in 

amorphous media if the formation length becomes large relative to the scat- 
tering mean free path of the electron. Nowadays, these effects are confirmed 
in experiment both for crystals and for amorphous media and are widely 
discussed. 

Several theoretical approaches extending the classical treatments of Refs 
[0, |2| have been proposed. Migdal [[| has used the Focker-Planck equation 
and quantified the results for amorhous media so that the corresponding 
effect is now known as Landau-Pomeranchuk- Migdal (LPM) effect. The 
polarization of the medium suppresses the soft photon emission as well (Ter- 
Mikaelian effect, see [|]]). The Kharkov group || has applied the path- integral 
technique to treat both effects simultaneously. In connection with recent 
SLAC experiments || on LPM effect, Blankenbecler and Drell have pro- 
posed to use at high energies the solution of the Klein-Gordon equation with 
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account of higher order corrections in expansion of the phase of the wave 
function in inverse powers of the initial momentum. The purely diagram- 
matic approach has been advocated in the series of papers done in Orsay 
[HI which extended the treatment of LPM effect to the non-Abelian case as 
well, in search for analogous effects in QCD started by Gyulassy and Wang 
|J. Somewhat special treatment using the Schrodinger equation and con- 
sidering compact quark- ant iquark-gluon systems has been proposed by B. 
Zakharov [[UJ who criticized some approximations used in papers ||. Soon 
it was shown [[□]] that with account of some additional terms omitted in |§ 
the two approaches are equivalent. Recently, Novosibirsk group [E| carefully 
considered different limiting cases of LPM effect in QED confronting them 
to experiment ||. Even though being common in spirit and close in final 
results, these approaches use different technique and different models of a 
medium so that sometimes it is hard to judge the correspondence between 
them. At the same time, some limiting cases are preferable to treat by either 
one or another method. 

Here, we would like to fill in one link in this chain of proposals and to 
show that for the Abelian case the solutions of the Klein-Gordon equation 
obtained in Ref. directly correspond to results of summing up the series 
of Feynman graphs considered in Ref. [HJ in the high energy (p — > oo) limit 
i.e. up to the terms of the order of 0{p" 1 ). In doing this, we show that the 
problem can be stated in terms of the post-eikonal approximation because 
of factorization of sums of Feynman diagrams in this limit. To simplify 
the formulas, we consider charged scalar particles since spin effects can be 
incorporated in a straightforward manner (see J?j). 



2 High-Energy Wave Function in the Post- 
Eikonal Approximation 

Let = {p°,p) = {\/p 2 + m ' 2 i 0)P) be the incoming four-momentum of a 
charged scalar particle, and 

(j)( x ) = e ip - x J d 3 qe^ s 4>(p; q) (1) 
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Figure 1: A Feynman tree diagram for the emission of n photons from an 
energetic particle of momentum p 

its energy eigenfunction with energy p° in the presence of a static source 

A°(x) = V{x) = f d 3 ke il:S v{k), (2) 



which transfers a total amount of momentum q to the particle. We assume 
the momentum k = (k, k z ) provided by the source at each interaction to be 
much less than the incoming momentum p (more exactly, k 2 <C 2p ■ k), so 
that the post-eikonal approximation 



1 1 1 / 1 k 2 1 



(p + k) 2 -m 2 + ie 2p-k + k 2 -ie 2p\k z -ie 2p (k z - ie) 2 
can be applied to the particle propagators of the Feynman diagrams. The first 



(3) 



term in the expansion constitutes the familiar eikonal approximation [fT3[] , but 
the second term is also needed to describe the Landau-Pomeranchuk-Migdal 
(LPM) effect. Our aim is to calculate (j)(x) in perturbation theory to an 
accuracy of order of 1/p in all orders. 

The triple interaction vertex is simply 2pv(k) in momentum space, and 
the seagull vertex |w (kijvfa) is down by a factor of p compared to the triple 
interaction. The amplitude in the tree approximation (see Fig. 1) is given by 



(f)(p;q) 


n=0 


4>n(p; q) 


= J (fl<l 3 kv(k,)jp n (p;h 


Pnip'i k±, ■ ■ ■ , k n ) 


= n 1 2p 

t=i (P + Ki) 2 - m 2 + ie 




i 



(4) 
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where the ellipses represent seagull contributions, and 4>o(f, 0) = 1- Note that 
the repeated appearance of v (fc)'s in (Q) is a result of multiple interactions 
with the same source v; it does not necessarily imply the presence of n distinct 
scatterers, though that can be accommodated. Note also that P(p; ki, ■ ■ ■ , k n ) 
is not symmetric in its variables k^. 

The perturbative expression for the wave function (|]) is 



x 



£e*> n (x), 



n=0 



M*) = / II e ih - s v(h)) Pn(p; h, ■ ■ ■ , k n ). (5) 
J 1=1 

In the eikonal approximation where the 0(k 2 /2p) terms in (|3]) are ne- 
glected, P n becomes 

n i 

p(p)( P - kl ,...,k n ) = (-ri[——. (6) 

If the correction terms in (H) are included to compute P n = P^ + Pjp/p + 
0(1 /p 2 ), then the first post-eikonal contribution is 

1 n f>2 

The corresponding contributions to n will be denoted by (p^ and (j>^ re- 
spectively. 

We shall use (f(ki,k 2 ,---,k n )) to represent the permutation average 
of any function /. Thus, for example, (f(ki,k 2 ,k 3 )) = (f(ki, k 2 , k 3 ) + 
f(h, h, k 2 ) + f{k 2 , h, k 3 ) + f(k 2 , k 3 , h) +f(k 3 , h, k 2 ) +f{k 3 , k 2 , h))/3\. We 
may and shall replace P n (k\, ■ ■ ■ , k n ) in (|5]) by its symmetric form {P n (ki, • • • , k n )), 
because this allows us to use the eikonal factorization formula (see also 
Appendix A) 

to compute the eikonal wave function: 

0i O) O?) = - } (-iXaW) n , 
nl 



4 



Xo(x) = Xo(z,b) = -i \ d i ke lkx v(k) 
(o) (x) 



k 



ie 



V(z',b)dz', 



cxp 



ip x + ipz — ixo{x) 



(9) 



When the post-eikonal terms in (3) are included, it is no longer clear how 
factorization and summation can be carried out. In order to get an idea how 
the post-eikonal contributions can be organized to yield factorization, let us 
look at the second-order contribution to cj)^ l \x): 



= f d%d%e i ^> s v(kMh)(P 2 [1 \p;h,k 2 )), 



(P 2 (1) (p; h,k 2 )) 



+ 



4 V (ku - ie) 2 k lz + k 2z - ie k lz - ie (k lz + k 2z - ie) 2 



+ 



k 2 1 1 (h + k 2 ) 2 



(k 2z - ie) 2 ku + hz 
1 1 



+ 



k 2z - ie (k u + hz ~ie) 2 



2 k lz + k 2z 



ie 



(10) 



The last term comes from the seagull diagram and the first four terms come 
from the post-eikonal contribution of the n! = 2! diagrams of the type shown 
in Fig. 1. The coefficient for k\ is 



1 



1 



4k lz + k 2z -ie\(k lz -ie) 2 
1 1 1 



+ 



ku 



+ 



-iz ~ 



1 



hz + k 2z - ie 



4 (k lz - ie) 2 k 2z - ie 



(11) 



Similarly, the k\ coefficient is —ik 2z — ie)~ 2 {k\ z — ie) -1 /^, and the ki ■ k 2 
coefficient is — [(ku — ie)(k 2z — ie)(ki z + k 2z — ie)]~ 1 /2. Combining this last 
coefficient with the seagull term in (|T0D reduces k\ ■ k 2 to the transverse dot 
product ki • k 2 . In short, the rather complicated expression in ( ]TD| ) can be 
reduced to the vastly simpler form 



[P?(v\ h,k 2 )) 



1 1 1 



1 



+ 2 



2! 2 ku - ie k 2z 
ki • k 2 

ku + k 2z - ie 



k\ 



ie \ k\ z — ie k 2z — ie 



(12) 
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Moreover, it is shown in Appendix A that this simplification occurs at all n, 
so that 



<P«(p; fci, •••,*»)> 



1 i ™ i 

TT 1 



i=l 



E 



it? 



+ 2E 
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(13) 



With this formula, the post-eikonal component of the wave function can be 
factorized into 



;r 



z,b) 



(n-ll! 



X2(^,b) 
Xi(^,b) 



- / d 3 ke ik - £ v(k) 



1 

2i- c 



v -*Xo 



i n-2/ 



(A; 2 -ze)2 
- / d 3 kd 3 k'e i{U ^ £ v(k)v(k' 



(n-2)! 

^'V 2 Xo(^',b), 
kk' 



-*XiJ» 



^' (V ± Xo(^,b)) 5 



(k z - ie)(k' z - ie)(k z + k' z - ie) 

(14) 



This factorization allows the sum over rz to be carried out to yield, to accuracy 

Ob- 1 ), 



cxp 



-ip°x + ipz - ixo(%) - i- (Xi(^) + %(^)) 

P 



(15) 



This expression agrees with the wave function obtained by Blankenbecler and 
Drell [[F] by solving the Klein-Gordon equation 



(E - V) 2 + V 2 - m 2 



0. 



(16) 



They looked for a solution of the form 4>(x) = exp[z$(a;)] accurate to order 
l/p in $(x), and found ( |I5] ) to be the solution. As shown in Ref. M the 
extension to spinor particles is straightforward. 

In summary, we have demonstrated that factorization of sums of Feynman 
diagrams does occur even in the post-eikonal approximation, in such a way 
to enable the perturbation series for the wave function to sum up to an 
exponential form, a form that agrees with the one obtained directly by solving 
the Klein-Gordon equation to accuracy 0(p _1 ) in the phase of the wave 
function. 
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3 Outgoing Wave Function 

In the previous section we have computed the energy eigenfunction <f>(x) 
with incoming momentum p^. In a similar way we can compute the energy 
eigenfunction (j)'{x) with outgoing momentum p ,fl : 



ip-x I d 3 q ' e i<! -asp (p>iF), 



& (p'-q 1 ) 



n=0 



J (f[ d?k< v(k'^ P< n {pi- k[,-.., k' n )8 (± k[ - q'^j 



2p' 



\ (p' — K-) 2 — m 2 + ie 



+ ■ 



(17) 



j'-i 



Since the transverse component p' of p' is not necessarily zero, the post- 
eikonal expansion fl3|) must be modified to read 



(p' - k') 2 - m 2 + ie 2p'-k'-k' + 



ie 



1 k' -2p'-k' 1 

+ 



2p' \ A£ + ie ' 2p' (k' z + ie) 2 
Following very similar arguments as before, we finally obtain 



(18) 



exp 



-ip' ■ x - ix'o(x) - i— (x'i(x) + 1X2^)) 



Xa(x) = i / d 3 ke lkx v{k) 



X' 2 (^ b ) 
X'i(^b) 



d 3 ke ik - s v(k) 



k z + ie 
k 2 



V(z',b)dz' 

1 r 00 



(k z + ie) 2 2 
d?kd 3 k'e l{m ' y *v(k)v{k') 



dzV 2 X ' (z',b), 
kk' 



(k z + ie)(k' z + ie)(k z + k' z + ie) 
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+ - /" d 3 ke^ s v(k) 



(k z + iey 

1 Z" 00 r o n 

= 2J Z dZ [( V ^o(^b)) 2 + 2p'-VxX / (/,b)]. (19) 

Let us note that the term linear in v appears in x'i because the transverse 
momentum p' differs from zero while Xi contains quadratic in v terms only. 



4 Multiple Scattering and Bremsstrahlung 

The on-shell scattering matrix element, with incoming momentum p^ and 
outgoing momentum p" 1 , is given by 



m fi = I d 4 x<f)'*(x)V(x)e ip - x = J d 4 xe~ ip - x V(x)(f)(x) 

= E / IKM&) )PL{p'\K---,k n ){2^[ P '-p-Y,h) 

n=l J \i=l ) \ i=l / 

oo„/n \ / n\ 

= E / IT rf3 M^) Pn(p;fcl,---,fcn)(27T) 4 5 4 p'"P-E fc i > 
n=l J \i=l ) \ i=l / 

(20) 

with = because the source is static. 

If, as a result of the scattering, a photon of four-momentum r M and po- 
larization vector £*(r) is emitted from the scalar particle, then the matrix 
element is 



Mfi(r) = —ie J d A x(f)'*(x)e*(r) ■ d 4>{x) 

oo / n \ 

= E <? V*(r) • 2p + 2 £ - r Mp; q) 

m,n=0 \ i=l I 



(m+n \ 
p' + r-p- h) ■ (21) 



The results obtained in the previous sections for the wave functions enable 
us to claim that up to accuracy 0{p~ 1 ) ) the treatment of Abelian LPM ef- 
fects by Blankenbecler and Drell who consider the Klein-Gordon equation 
directly, corresponds to that of R. Baier et al pfl who use the diagrammatic 
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approach. At high energies and for finite targets the first approach can be- 
come preferable because of the usage of simplified expanded propagators and 
direct treatment of spatial evolution. 

It could become especially simple in case of thin targets, where emission 
at single scattering with some corrections due to double scattering (see Ap- 
pendix B) dominates, and formulas for P± and as given by flT3"| ) are 
exploited. The slightly corrected Bethe-Heitler regime is at work. 

For finite targets, there are three lengths important in the problem. Those 
are the target thickness I, the mean free path l m and the formation length If. 
For soft photons, If ~ 2 r y 2 ju). Here 7 = p°/m, u = r° is the photon energy. 
For screened Coulomb fields in QED, Z" 1 = Anr 2 Z 2 \n(183/Z 1 ^ 3 ), where n is 
the density of the scattering centers, r e = a/m « 2.8 fm. Let us consider 
the case of very thin targets and high energy electrons when I <C l m <C If. 
We show that the decline from the Bethe-Heitler formula is determined by 
the target thickness in units of the mean free path. 

Following Ref. 0, one can express the total intensity of the radiation 
I(tot) as the Bethe-Heitler intensity I(BH) m AT (where T = irl/3l m ) sup- 
pressed by a form factor F: 

I{tot) = I(BH)F. (22) 

From formulas (9.16) of Ref. [0, it is easy to find out that for thin targets in 
the soft photon limit the form factor is given by 

- 1.5T. (23) 

The correction is small for T small enough and vanishes linearly with target 
thickness. 

At first sight, it seems that the small sizes of hadronic targets favor this 
limiting case for non-Abelian effects. However, this statement requires fur- 
ther study since the nuclear target thickness, in units of the mean free path 
T nuc i (which is analogous to T), could be rather large. Then, inspired by eq. 
( |2"3"D of QED, one would guess that the QCD suppression can become strong, 
i.e., nuclear LPM effect essential even though the target thickness in units of 
the formation length is very small. 
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5 Conclusions 



Thus we have shown the equivalence of two approaches to the problems 
of scattering and radiation of high-energy electrons traversing the dense 
medium. One of them |7| deals with the solution of the Klein-Gordon equa- 
tion, and another one || with sums of Feynman graphs. Even though the 
initial formulas have quite different forms, they lead to essentially the same 
expressions for the solutions of the above problems. This is due to the fact 
that the sum of Feynman diagrams calculated at 0{p^ 1 ) accuracy factorizes 
still in such a way that it can be represented by the corresponding post- 
eikonal expansion of the phases of the wave functions. 

Besides, we have argued that the nuclear LPM effect can be essential even 
though the size of nuclear targets in metric units is very small. 
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A Factorization in the Post-Eikonal Approx- 
imation 

The aim of this appendix is to prove eq. (|13|). It would be useful first to 
review how (§) is arrived at. 
The integral 

gives the left-hand side of (|j) before taking the permutation average, where 
the hyper-triangular integration region R(12 ■ ■ ■ n) is defined to be {0 > t\ > 
t 2 > ■ ■ ■ > t n > — oo}. The permutation sum on the left of (|§P can be 
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obtained by summing the integral over all permuted regions, whose union 
is the hyper-rectangular region {0 > tj > — oo}. Upon integration over the 
hyper-rectangle and a division by n\, one obtains the right-hand side of (^) 
and hence the eikonal formula. 

The k 2 m coefficient of P^Kh, ■■■,k n ) can be seen from (|^) and (|7|) to be 

idL^ih,---,^. (25) 

Moreover, this relation persists upon permutation averaging. Thus the k 2 ^ 
coefficient of (P^) can be obtained by differentiating both sides of @ with 
respect to k mz . The result is the one given in fll3|) . 

It also follows from (§) and (0) that the 2k^k m coefficient of P^'(ki, ■ ■ ■ , k n ) 
can be obtained by differentiating P^ with respect to k pz , where p is the 
larger of the two numbers i and m, viz., the one that stands to the right 
of the identity permutation (12 • • -n). Upon permutation, this relation still 
holds provided p is taken to be the number form that stands to the right 
of the other. The permutation sum will now be divided into two sums, one 
with £ standing to the right of m, and the other with m standing to the 
right of I. In the first case, the total integration region is a product of the 
triangular region {0 > t>m >t^> — oo} with the hyper-rectangular regions 
{0 > ti > — oo} of the remaining n — 2 variables. The result is 

n 1 \ 1 1 

n — — — — — ( 26 ) 

m k% z l€ J k mz zc k mz -)- kg z Z6 

Upon differentiation with respect to k pz = ki z it produces an additional factor 
of -(ktz + k mz - ze) _1 . 

In the second case, m and I are interchanged, but the differentiation with 
respect to k pz = k mz still produces the same additional factor as before. 
Adding up these two cases, we obtain ([]) multiplied the additional factor 
— kn z + k mz — ze) -1 , which is then equal to the coefficient of 2kg ■ k m in flT3p. 
As before, when combined with the seagull contributions, 2ke ■ k m is reduced 
to 2k^ • k m as shown in 
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B Thin Targets 



For thin targets, the photon emission at single scattering (Bethe-Heitler 
regime) becomes prevailing, though there are corrections due to double- 
scattering processes. We write down some formulas for this case. 

According to eqs (|14D , (|l5l), (|T9|) , (|2~0"D , and fl2l]) , the matrix element of 
soft photon emission by a scalar electron, scattered once with longitudinal 
momentum transfer much less than p, is given by 



M l oc -i ( d 3 xs*(r) ■ plxo + ~{x\ + ixl)} 
J p 



(27) 



where Xj = Xj + x'j-, an d only terms linear in v are kept in xi an d x'\- It 1 



is 



easy to get 



x\ + ixl 



xi 



2tt J rf 2 kw(k,0)e ik ' b 



- / d 3 kv{k)e 



ik-x 



+ 



(28) 
(29) 



These expressions follow directly both from Feynman diagrams and from the 
solution of the Klein-Gordon equation obtained in Ref. 0. 

The matrix element of photon emission at double scattering is given by 
eq. (21) for m + n = 2. The two terms with m = and n = describe 
emission before or after the scattering, while the term with m = n = 1 
corresponds to photons emitted between two scatterings. It is easy to check 
that 1/p-contribution contains the following factor in the integrand: 



1^ oc k-k' 



1 



+ 



- ni6(k z )k' 2 

- m5{k' z )k 2 

- ni5(k 



(k z - ie){k' z - ie)(k z + k' z - ie] 
1 

{k g + ie){k' z + ie)(k z + k' z + ie) 
1 1 



{k' z -ief 
1 

p'k' 



(K+iey 



ni5(k' z 



i 

Jk~+W 
p'k 



(k z + ief 



(30) 



12 



This expression can be also obtained both directly from propagators in Feyn- 
man graphs and from phases of the wave functions. Its structure is clear from 
above formulas. 
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